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ABSTRACT 

Recent observations of cosmic ray electrons from several instruments have 
revealed various degrees of deviation in the measured electron energy distribution 
from a simple power-law, in a form of an excess around 0.1 to 1 TeV energies. An 
even more prominent deviation and excess has been observed in the fraction of 
cosmic ray positrons around 10 and 100 GeV energies. These observations have 
received considerable attention and many theoretical models have been proposed 
to explain them. The models rely on either dark matter annihilation/decay or 
specific nearby astrophysical sources, and involve several additional assumptions 
regarding the dark matter distribution or particle acceleration. In this paper 
we show that the observed excesses in the electron spectrum may be easily re- 
produced without invoking any unusual sources other than the general diffuse 
Galactic components of cosmic rays. The model presented here assumes a power- 
law injection of electrons (and protons) by supernova remnants, and evaluates 
their expected energy spectrum based on a simple kinetic equation describing 
the propagation of charged particles in the interstellar medium. The primary 
physical effect involved is the Klein-Nishina suppression of the electron cooling 
rate around TeV energies. With a very reasonable choice of the model parameters 
characterizing the local interstellar medium, we can reproduce the most recent 
observations by Fermi and HESS experiments. Interestingly, in our model the 
injection spectral index of cosmic ray electrons becomes comparable to, or even 
equal to that of cosmic ray protons. The Klein-Nishina effect may also affect 
the propagation of the secondary pairs, and therefore modify the cosmic ray 
positron-to-electron ratio. We have explored this possibility by considering two 
mechanisms for production of pairs within the Galaxy. The first is due to 
the decay of 7r s produced by interaction of cosmic ray nuclei with ambient 
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protons. The second source discussed here is due to the annihilation of the 
diffuse Galactic 7-rays on the stellar photon field. We find that high positron 
fraction increasing with energy, as claimed by the PAMELA experiment, cannot 
be explained in our model with the conservative set of the model parameters. We 
are able, however, to reproduce the PAMELA (as well as Fermi and HESS) results 
assuming high values of the starlight and interstellar gas densities, which would 
be more appropriate for vicinities of supernova remnants. A possible solution to 
this problem may be that cosmic rays undergo most of their interactions near 
their sources due to the efficient trapping in the far upstream of supernova shocks 
by self-generated, cosmic ray-driven turbulence. 

Subject headings: cosmic rays — Galaxy: general — ISM: general 



Introduction 



Measurements of the energy spectra of cosmic ray (CR) species are of a great impor- 
tance for understanding the physics of Galactic CR sources (such as pulsars or supernova 
remnants), as well as for constraining the internal structure of Milky Way, since this struc- 
ture (topology and intensity of the Galactic magnetic field, profiles and energy distribution 
of different Galactic photon fields, distribution of interstellar gas and dust, etc.) determines 
the spatial and energy evolution of the injected charged particles which propagate through 
the interst ellar medium. In addition, as discussed in several papers by a number of au- 
thors (e.g.. lJungman et al.lll996t ICheng et al.ll2002l ). annihilation or decay of a dark matter 
(Kaluza-Klein particles or supersymmetric WIMPs) may also imprint some signatures in the 
observed spectra of CR electrons within the GeV-TeV energy range. Conseque ntly, the most 



recent observations by the ATIC, PAMELA, Fermi, and HESS experim ents ( jChang et al. 



20081 : lAdriani et alJl2009t lAbdo et alJl2009t lAharonian et al.ll2008l . l2009al . respectively) have 
generated significant interest on this topic in astrophysic and particle physics communities. 

The most sophisticated framework for analyzing the CR propagation within the G alaxy 
is provided by the GALPR OP model presented first by iMoskalenko fc Strong! (119981 ) and 
Strong fc Moskalenkd (119981 ). This model assumes injection of a power- law electrons and 
nuclei, and follows their spatial and energy evolution under the influence of different radia- 
tive processes (Coulomb losses, synchrotron and IC cooling, proton-proton collisions, etc.), 
taking also into account the relevant interactions of charged particles with the interstellar 
turbulent magnetic field (with the assumed Kolmogorov spectrum) in a quasi-linear approxi- 
mation regime. The model can successfully explain many find ings regarding the hadronic CR 
spectrum and its composition (IMoskalenko et al.ll2002u2003l ). as well as the observed Galac- 
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tic di ffuse 7-ray emission (IMoskalenko fc Strongll2000l ; IStrong et al.ll2000l . I2004J : iPorter et al. 
20081 ). However, the model predicts also the decrease of the CR positron-to-electron ratio 
with particles' energy in the GeV-TeV range, in a disagreement with the observational in- 
dications. This prediction was made under the working hypothesis that bulk of the Galactic 
e ± pairs are created in the collisions of relativistic CR proton s with ambient gas, and th e 



subsequent decays of the generated pions (see the discussion in IMoskalenko fc Stronglll998l ). 
One should note that the previous (prior to 2008) measurements regarding this issue were 
restricted to electron energies E e < 10 GeV, and as such could be seriously affected by the 
charge dependence of solar modulation. On the other hand, the most recent PAMELA ob- 
servations, reachin g E P ~ 100 GeV en ergies, confirmed the increasing positron fraction in 
the CR spectrum (lAdriani et al.ll2009l ). 



The other challenge to the 'standard' model of CR propagation came from the observa- 
tions of ATIC collaboration, which reported a sharp pile-up around E e ~ 0.5 TeV above the 
power-law sp ectrum both obser ved at lower energies but also emerging from the GALPROP 
calculations (IChang et al.ll2008l ). Reality of this sharp spec t ral feature has been que stioned 
by the Fermi and HESS experiments (see lAbdo et al.ll2009l ; lAharonian et al.ll2009al ). which 
show a much smaller and broader excess over the best fit power-law continuum J e (E e ) oc E~ 3 . 
In this context one has to keep in mind that due to rapid radiative losses of the TeV-energy 
electrons, their spectrum measured in the Solar System may be possibly domi nated by a 



few local sources (most likely nearby pulsars such as Vela or Geming a; see, e.g.. IShenlll970 



Nishimura et al.l Il980l ; lAharonian et al.l Il995l ; iKobayashi et al.l 12004 ) , and therefore may be 
more complex than a featureless power-law continuu m. In particular, it may reflect the 



non- s tationary and stochastic nature of such sources (IPohl fc Espositdll998l ; iGrasso et al. 



20091 ). More interestingly, the substructure around 0.5 TeV observed by ATIC was argued 



to be con sistent with that e xpected from the annihilation of the Kaluza-Klein dark matter 



particles (IChang et al.ll2008l ). Both 'local pulsar' and dark matter scenarios were claimed to 



Pohl 


2009; 


Grasso et al. 


2009) 



Obviously, the dark matter interpretation of the ATIC, Fermi, HESS and PAMELA 
results is of a great interest. However, in order to fit the collected data in the framework of 
this model, large though arbitrary 'boost factors' have to be invoked for the dark matter an- 
nihilation/decay fluxes, which are only roughly justified by a possible non-unifo r m (clumpy) 



distribution of the dark ma tter near the Solar Syste m fe.g. IChang et al.ll2008l; lElahi et al. 
2009 ; Hooper fc Zurek 2009 ). In addition, as shown by Profumo fc Jeltema (j2009l ). the dark 



matter scenario for the positron excess conflicts with the observ ations of the extragalac - 
tic background radiation in the X-ray/7-ray energy range (see also lBelikov fc Hooperll2009l ). 
Similarly, the 'local pulsar' interpretation also relies on several model assumptions (regarding 
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particle acceleration in relativistic outflows), which are not observationally verified yet. 

In this paper we explore the possibility of explaining the aforementioned data sets by 
a simple model for the generation and propagation of CRs in the Galaxy (and the vicinity 
of Earth), without invoking any new/unconventional sources of the TeV-energy electrons 
or positrons. This scenario is not intended to compete with the complexity of the GAL- 
PROP model. It is intended instead to point out several effects which, even though being 
standard and at some level inevitable, might have been underestimated or overlooked in 
the previous analysis. The new physics involved here is that we emphsize the importance 
of the Klein-Nishina (KN) suppression of the inverse Compton (IC) scattering cross section 
for ultrarelativistic electrons and positrons interacting with the Galactic starlight and other 
photon fields. We evaluate the modification of the primary electron energy spectrum which 
results from this new aspect, as well as that of secondary pairs arising from the annihilation 
of the Galactic 7-rays on starlight and from CR proton interactions with the interstellar 
medium (ISM) . In § 2 we discuss our model and present the kinetic equation describing the 
propagation of CR electrons (as well as protons and positrons). We first describe the KN 
effect qualitatively and then compare the CR electron spectra obtained from detailed solu- 
tion of the kinetic equation with observations. In § 3 we consider various mechanisms for 
production of secondary e ± pairs and derive their spectra using the same propagation model, 
comparing again the resulting spectra with those observed by PAMELA. A brief discussion 
and summary is presented in § 4. 

2. Primary CR Electrons 

We start this section with the description of the relevant interactions of CR electrons 
(accelerated and injected into the Galaxy by, presumably, supernova remnants) with the ISM 
photons (mainly starlight), magnetic field, and turbulence. We then describe the resultant 
spectra of these 'primary' electrons and compare them with the observations reported by 
ATIC, Fermi, HESS, and other experiments. 

2.1. Interactions of CR Electrons 

The propagation of relativistic electrons injected into the ISM is determined by two basic 
interactions: radiative cooling and interactions with plasma turbulence. The latter causes 
diffusion in space (determining the rate of the escape of electrons from the Galaxy) and 
diffusion in energy (determining the rate of the acceleration). For the electrons with energies 
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above GeV, the radiative cooling is mainly via the IC scattering on ambient photon fields and 
due to the synchrotron emission in the Galactic magnetic field. The relevant photon fields 
are the Cosmic Microwave Background (CMB) radiation with the energy density u cm b ~ 
0.26 eV cm -3 , the Galactic starlight, and far infrared photons from the dust emission. The 
latter two are expected to dominate over the CMB (i.e., u sta r, Udust > Ucmb) in the inner 
part s of the Galactic d i sk (r 5, 10 kpc from the center) by a factor of at least a few (see, 



Strong et allboool : [porter et alJkooii iMoskalenko et alJliooj : IPorter et allbood ). For 



e.g. 

an electron with energy E e , the characteristic cooling timescale in the Thomson (T) regime 
is therefore given as 



Trad, T 



3m 2 c 3 
4cr T Utot E e ' 



with 



B 2 



u t ot = u cm b (1 + + ^— and £ = 



57T 



Udust U s tar 
Ucmb 



From this one gets r rarf T ~ 10 9 (Ee/GeV)' 1 (1 + 0.1 B 2 G + ^)~ 1 yr, which, for the illustrative 
(though expected) Galactic magnetic field B^q = B/fiG ~ 3 and £ ~ 10, leads to r m d,r — 
100 (i? e /GeV) _1 Myr. This estimate breaks down at high particle (and/or photon) energies 
where the KN effect reduces the IC cross section. This happens when the target photon 
energy in the electron rest frame exceeds electro n rest mass, which, in the o bserver frame, 



translates to e > m 2 e c 4 /AE e ~ 65 (EjGeV)- 1 eV flBlumenthal fc Gouldlll970h . 



For the evaluation of the interactions of CR electrons with magnetic turbulence we use 



the quasi-linear approximation for the particle- wave interactions (see, e.g.. ISchlickeiserl 12002 



and references therein) and assume a Kolmogorov spectrum of turbulence as a superposition 
of magnetohydrodynamical (MHD) waves. The corresponding spatial diffusion timescale 
may be estimated as r esc — 3£ 2 /cX^E e ), where i is the linear scale of the system and the 
particle mean free path X(E e ) ~ ( r^ 3 Xmax- Here £ = (B/5B) 2 is the ratio of energy densities 
stored in the large-scale ('unperturbed') and turbulent magnetic fields, r g = E e /eB ~ 3 x 
10 12 (i? e /GeV) B~q cm is the electron gyroradius, and X m ax is the maximum wavelength of 
the turbulent modes. This gives r esc ~ 10 7 C" 1 B l J 3 (£/kpc) 2 (A ma:E /kpc)- 2 / 3 {Ec/GeV)- 1 / 3 yr 
which, for the ISM parameters B^ G ~ 3, £ ~ 1, and X ma x — lkpc, simplifies further to 
T esc — 10 (£/kpc) 2 (i? e /GeV) _1 / 3 Myr. In other words, an electron with energy E e travels the 
distance I ~ 3 (Ee/GeV)^ 1 ^ 3 kpc within ISM before loosing its energy via radiative cooling 
(for £ ~ 10). Hence, high-energy CR electrons (E e > 10 GeV) detected near the Earth are 
suppose d to origina te from local (£ < 3 kpc) region and recently operating (t < 100 Myr) 
sources flShenlll970l ). 



It should be noted in this context that several observational findings (regarding, e.g., 
secondary-to-primary ratios of some CR elements) are better interpreted in terms of the 
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particle diffusion shaped by the ISM turbulence characterized by th e Kraichnan energy 



spectrum, rather than of the Kolmogorov form anticipated above (see iPtuskin et al.l 12006 



Strong et al.ll2007l . and references therein). If this is the case indeed, then X(E e ) ~ £ r x J 2 xUax 
and r esc ~ 500 (£/kpc) 2 (E'g/GeV)" 1 / 2 Myr, and hence the conclusion regarding the local 
origin of > lOGeV energy CR electrons is even strengthened. 

On the other hand, the order-of-magnitude estimate presented above should be taken 
with caution, because of the uncertainty in the value of X max , which is not a directly mea- 
sured quantity, but is only expected to be roughly of the scale of the Galactic disk thick- 
ness, ~ 1 kpc. In addition, the use of the Kolmogorov/Kraichnan spectru m for ISM turbu- 



lence is only partl y justified by observations/theoretical models (see, e.g.. ICho et al.l 12003 



Strong et al.l 120071 . and references therein). Finally, the anticipated quasi-linear approxima- 
tion is appropriate for modeling of the particle diffusion along the magnetic field lines (and 
thus rather within the Galactic disk), and is not expected to describe properly the prop- 
agation of particles accross the magnetic field lines (i.e., the escape of particles from the 
Galactic disk). 

The spatial diffusion of CR electrons is accompanied by their diffusion in the momen- 
tum space, leading to a net acceleration of particles on the characteristic timescale r acc ~ 
3 X(E e ) I ' c0j„ where xr . is the vel ocity of turbulent MHD modes in the units of speed of light 
(e.g., iBlandford fc Eichlerl 119871 ). For a low-beta plasma, the MHD wave velocity is equal 
to the Alfven velocity c(3 sc ~va — -B/(47rm p n ism ) 1 / 2 ~ 2 x 10 5 B^ G (n ism /cm" 3 ) -1 / 2 cms -1 , 
where ni sm is the number density of the ambient plasma. Thus, for a Kolmogorov spec- 
trum of the ISM turbulence r acc ~ 10 11 ( B^ /3 (A max /kpc) 2 / 3 (n lsra /cm- 3 ) {EjGeV) 1 / 3 yr. 
Assuming again values B^q ~ 3, ( ~ 1, X max — lkpc, and n ism ~ 1cm -3 , this simplifies fur- 
ther to r acc ~ 10 4 (i? e /GeV) 1//3 Myr, which indicates that ultrarelativistic electrons undergo 
little turbule nt acceleration with in ISM before they cool radiatively or escape (see in this 
context, e.g.. IPtuskin et al.ll2006l ). Thus, radiative cooling and diffusive escape are the main 
processes controlling evolution of CR electrons. Again, in the case of the Kraichnan form of 
the magnetic turbulence, the same ISM parameters give r acc ~ 400 (-E e /GeV) 1//2 Myr, which 
is still much longer than the radiative cooling timescale of electrons with E e > 10 GeV. 



2.2. Transport Equation 

In what follows we assume that CR electrons are injected by numerous sources through- 
out the ISM which has a smooth and slowly varying distribution of gas, photons, and mag- 
netic fields relative to the relevant interaction scales discussed above. As indicated by the 
small mean free path calculated above, the electrons undergo multiple scattering before any 
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other interactions so that they acquire an isotropic pitch angle distribution. As a result, for 
the relevant scales of about few kpc in the vicinity of the Earth we can use the homoge- 
neous and isotropic approximation in describing the transport of the CR electrons, where the 
spatial diffusion can be represented by an overall (energy dependent) escape term. Hence, 
ignoring turbulent acceleration, the propagation of the GeV/TeV CR electrons within the 
local ISM can be described by the following kinetic equation: 



dn e {E e 
dt 







dE e 



E e njE e 

T~loss 



nlE e ) 



+ QlE e ) 



(2) 



(see, e.g.. IPetrosian &: Liull2004j . and references therein), where Q e (E e ) denotes the injection 
rate of electrons, and Ti oss the total energy losses timescale. The steady-state solution to 
this equation reads as 



nlE e 



Tioss i,E e 

E P 



dE' e Q e (K) exp 



l - dE'l r loss {E'l) 



(3) 



For a broad, e.g. a power-law-type injection function, a very rough but illustrative approx- 
imations to this solutions are n e (E e < E ) ~ r esc x Q e and n e (E e > E ) ~ Ti oss x Q e , 
where at E we have Ti oss (E ) = T esc (E ). Note that for T[ oss ~ r rarf T and r esc as spec- 
ified above (£ ~ 10, and i ~ 3 kpc, roughly the vertical scale of the Galactic disk) one 
has E ~ 1 GeV. Thus, from the precisely known transport timescales and the observed 
electron flux J e (E e > GeV) oc E~ 3 one can set direct constraints for the electron injection 
function Q e (E e ). We note in this context that non-relativistic shock waves associated with 
Galactic supernova remnants (SNR s) are expected to be the pr imary source of the observed 
high-energy (> 10 GeV) CRs (e.g.. iBlandford fc Ostrikerlll978l ). This anticipation seems to 
be confir med by the most recent X-ray and 7-ray observations regarding several remnants 



(see, e.g., iVinki 120081 . and references therein). However, the injection spectrum of CR par- 
ticles (both leptons and hadrons) may be substantially different from the energy spectra 
of freshly accelerated electrons and ions at SNR shocks, due to a complex convolution of 
particle cooling, transp ort, and CR-driven magnetic amplification pro cesses in vicinity of 
non- linear shock waves (IBlandford &: Ei chlerl 1 1 9 8 71 ; ICaprioli et al.ll2009l ). 



In this paper we restrict the analysis to the CR electrons not affected by the solar 
modulation, i.e. the ones with energies E e > 10 GeV, for which the escape effects can be 
neglected. Therefore we omit the escape term in the equation [2], and obtain 



n e (E e ) 



Tic 



E P 



dE' e QIK 



(4) 



so that for the injected function Q e (E e ) oc E e Se and ti 
electron spectrum will also be a power- law, n e (E e ) cx E~ St 



- t oc E e 1 the observed 
Hence s e = 2 will be required 
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by the observed electron flux J e (E e ) oc E~ z . However, if the dominant electron cooling is 
due to the IC scattering in the KN regime, the energy losses time scale goes roughly as 
t 1oss ~ T ra d,KN oc E^ 2 , which will give rise to n e (E e ) oc E~ Se+0 - 5 . In other words, within 
the energy range where the IC/KN losses are dominant, the steady-state electron energy 
distribution is expected to pile-up above that expected from e xtrapolation of the Thomson- 
regime spectrum, as discussed i n different contexts by, e.g., 



(119851): iDermer & Atovanl fl2002f ); Kusunose & Takaharal (120051 ) 



A 



laronian &: Ambartsumvan 
Moderski et all (l2005h . and 



Stawarz et al.l (j2006l ). These KN-related spectral pile-ups are more and more pronounced 
for flatter and flatter injection continuuno Interestingly, as already noted in the previous 
section, for the characteristic energy of the starlight photons e star — 1 eV (wavelengths 
X s t a r — lyum), this is expected to happen for E e > m 2 e c A /Ae sta r ~ 0.1 TeV, i.e. within the 



range of the claimed 'electron excess', assuming that the energy density of the starlight 
emission dominates all the other Galactic photon and magnetic fields (£ > 1). As a result, 
if the discussed KN effect plays a role, a relatively steep electron injection index s e > 2 is in 
fact required to account for the observed spectrum J e (E e ) oc E~ 3 . 



2.3. Energy Spectra of CR Electrons 



We now present a more rigorous treatment of the KN effect and show that it may be 
the primary cause of the observed deviation of CR electron spectrum from a simple power 
law. For this purpose we need a more detailed description of the Galactic photon fields and 
of the radiative cooling rate. We model the Galactic dust and starlight emission by the 
functions u(e) similar to the ones given by lPorter et al.l (120061 ). which peak for e ~ 0.015 eV 
and e ~ 0.5 — 2.0 eV photon energies with the maximum levels Ud us t and u star , respectively, 
such that the total (integrated over e) energy densities in these components are ~ 2 Udust and 
— 3 u s tar- The left panels of Figured] show variation with photon energy of the CMB, dust 
and starlight energy densities for several relative values of these densities. In each panel we 
show three curves: for three values of Ud us t and fixed u star = 0.3eVcm~ 3 (top); three values 
of u s tar and Udust = 0.3eVcm~ 3 (middle); three values of Ud us t = u s t ar (bottom). We also 
show the energy density of the magnetic field for B = 1 //G (solid horizontal lines), and 3 //G 
(dashed horizontal lines). 



Follwing iModerski et al.l (120051 ) , we approximate the radiative energy loss timescale for 



1 Note that the pile- up effects in the electron energy distribution resulting from the KN suppression are 
present at some level for any slope of the injection spectrum, unlike the analogous effects related to the 
synchrotron cooling alone (|Kardashevlll962n . which are present only for s e < 2. 
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ultrarelativistic leptons as 

Trad(Ee) ' 



2^3 



3m 
4a T E e 



8tt 



+ / deutole) f K N 



(±E e 



(5) 



where the total energy density of the Galactic photon fields (including CMB radiation) is 
J deutoti^)- Here the KN correction factor is taken to be of the form 



fKN{x) 



(1 + x)- 1 - 5 for x < 10 4 
far 2 (\nx-f) for x > 10 4 



(6) 



In addition, we consider the electron energy losses due to the Coulomb collisions and the 
bremsstrahlung process. The appropriate timescale of these can be approximated as, respec- 
tively, 



2E P 



3 InA m e c 3 <JT Ui S 



(7) 



where In A ~ 40 (jPetrosianl Il973l ) , and 

,(Ee) 



2tt 



3af s a T cn ism ^E e 



where otu — 1/ 137 is the fine structure constant and y^E e ) ~ In (2E e /m e c 2 ) — 1/3 ~ 10 
(jPetrosianl l200ll ) . We note that the above form of Tb rem includes electron- ion and electron- 
electron bremsstrahlung assuming completely unscreened limit with 10% fully ionized helium 
abundance. With such, the total energy losses timescale is 



r£(E e ) = r^ ul {E e ) + T^\jE e ) + r; a \{E e 



(9) 



The right panels of Figured] show the total energy losses timescales (multiplied by en- 
ergy) corresponding to the photon and B field energy densities the same as in the left panels, 
and TLi sm = 1 cm -3 . As evident, at low energies E e < 1 GeV the Coulomb and bremsstrahlung 
processes dominate electron cooling, since r coul ~ 50 (E e /GeV) Myr and r brem ~ 50Myr (for 
n ism — 1 cm -3 and }^E e ) = 10). At higher electron energies, the IC/T losses take over. 
However, for E e > 10 GeV the radiative cooling rates deviate from the ones characterizing 
the Thomson regime (which would be represented by horizontal lines on these plots) due 
to the KN effect. Note that these deviations are the strongest in the case of a large ratio 
Ustar/udust- This is because for dust emission in the far-infrared range, the dominant radia- 
tive cooling is still in the Thomson regime even for relatively energetic electrons. Therefore, 
the KN suppression for the optical target photons becomes important only for large values 
of the ratio u staT jud U st- This is exactly the reason why the KN-related features in the CR 
electron spectrum discussed here may remain unnoticed in the GALPROP calculations, even 



-10- 



though this code includes the exact prescription of the IC cross section, valid in both T and 
KN regimes. 

Using this radiative loss rate in the simplified version of the kinetic equation we 
obtain the energy flux spectrum J e (E e ) oc n e (E e ) of CR electrons. For the injection function 

Qe(E e ) we use 

Q e (E e ) = k e E~ S£ exp 



(10) 



with the normalization k e fixed so that [E^J e (E e )]E e =30GeV = 151.4 GeV 2 m -2 s _1 sr" as in- 
dicated by the Fermi data. In general, the model outlined above has seven free parameters, 
namely s e , E e ^ max , Udust, u star , B, n ism , and I. However, we fix for illustration t = 3kpc, 
Uism = lcm -3 , and E e ^ max = 2TeV, so that we are left with only four free parameters s e , 
Udust, u s tar, and B. We note in this context that the direct measurements of the Galactic pho- 
ton and magnetic field energy densities are difficult d ue to substantial foregrounds, and thus 



the a ssociated uncertainties are relatively large (e.g., ICrutcher et al.ll2003l ; lHauser & Dwek 



20011 ) . Below we explore the corresponding parameter space of the model. 



Figure[2] shows the energy spectra of primary electrons corresponding to two different 
injection spectral indices, s e = 2.0 and 2.2, and to the same choice of the values of the other 
three model parameters used in Figured! As evident, the expected KN pile- up effects are 
indeed present, being the most pronounced for flatter injection continuum, and for large 
values of the ratio u s tar/udust- The value for the ISM magnetic field have little effect on 
the results, as long as B < 10 /xG. One conclusion here is that different combinations of the 
parameters s e , u sta r, and Ud us t can lead to the observed electron spectrum J e (E e ) oc E~ 3 . 
However, our primary result is that it is relatively easy to account for a possible minor excess 
in the energy distribution of primary CR electrons over this power law in the 0.1 — 1 TeV 
energy range purely by the KN effect. 

In Figure[3] we compare the observed spectra from various experiments with one of our 
model calculations corresponding to a choice of model parameters appropriate for the average 
(local) ISM conditions, namely £ = 3kpc, n ism = lcm -3 , B = 3uG, Udust = O.leVcm -3 , 
Ustar = 3eVcm~ 3 , E e ^ max = 2.75 TeV, and s P , = 2.42 . The data points corr espond to d i fferen t 



measurements by ATIC (bla ck; Ichang et aljlgooj) PPB-BETS (yellow; iTorii et alihoosh. 



emul s ion cha mbers (magenta; iKobayashi et al.ll2004f). HESS (b lue and cyan; lAharonian et al. 
20081 . [2009a . respectively), and Fermi (red; lAbdo et all 120091 ). As evident, with reasonable 
parameters^ and the electron injection index the same as required for the Galactic CR protons, 



2 It is important to note in this context that the anticipated value of the starlight energy density, u star = 
3eVcm -3 , even though considered here as a 'reasonable' one, is still larger than that expected for the local 
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s e = Sp = 2.42 (see section 3.1 below), we can reproduce very well the latest, and most 
reliable observations by Fermi and HESS. It is also clear that the KN effect cannot account 
for the sharp feature claimed by ATIC observations. Encouraged by this simple and robust 
explanation for the primary CR electron spectrum, in the next section we also explore the 
influence of the KN effect on the expected spectra of secondary pairs with the goal of 
providing explanation of the rise with energy of the positron to electron ratio observed by 
PAMELA. 



3. Secondary and Tertiary Pairs 

In this section we address the question of the origin and spectrum of ultrarelativistic 
positrons present in the CR population, which are produced as secondaries in pair pro- 
duction processes. We consider three different sources of secondary e ± pairs and apply the 
same transport equation as above to determine their spectra in the ISM. From these we 
obtain the positron to electron ratio and compare it to the observation by PAMELA. 



3.1. Proton-Proton Pair Production 

The first source of secondary pairs we consider is due to the interactions of ultrarela- 
tivistic CR ions (primarily protons) with the ambient plasma. We assume that the Galactic 
sources of CRs in addition to electrons inject also ultrarelativistic protons at a constant 
rate Q P (E P ) oc E p Sp , which then propagate diffusively through the ISM and collide with 
cold protons. The appropriate timescale for the proton-proton interaction is roughly in- 
dependent of energy: r pp ~ (cn ^ m a pv )~ l ~ 30 (n^ m/cm -3 ) -1 Myr for the cross-section 



a 



3.4 x 10 



-26 cm 2 



ipp - 

[see, 



e.g., 



Kelner et al.l 120061 ). The diffusive escape timescale for 
CR protons is same as for electrons, namely r esc ~ 100 (Ep/GeV)^ 1 ^ 3 Myr (for £ ~ 3kpc, 
and rii sm ~ 1cm -3 ). This means that CR protons with E p > 30 GeV are in a slow cool- 
ing regime (i.e., we are dealing with a thin target case), so that the ISM proton energy 



spectrum can be approximated as n p {E p > 30 GeV) ~ r esc x Q P (E P ) oc E t 



-Sp-l/3 



(for the 



Kolmogorov turbulence; see the discussion in section 2.1 and below equation^]). Keeping in 



mind the observed CR proton flux J p (E p ) cx E p , the required injection spectral index 
should be then s p ~ 2.42. In addition, in this regime protons escape with most of their 



ISM (jstrong et al.ll200d [Porter et all [20061 : iMoskalenko et allbood: borter et al.ll2008h . As such, it should 
be considered as an illustrative model assumption, for which the analyzed KN effects are already of a major 
importance (see the related discussion in §3.4 and §4 further below). 
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energy and only a small fraction / of the carried flux goes into production of secondaries 
(e^ and neutrinos arising from 7r =1= decays) and 7-rays (from tt decay). In particular, one 
has / ~ r esc jr vv ~ 0.3 (.Ep/TeV)" 1 / 3 . Note that if the CR protons propagate through the 
ISM with some particular bulk velocity, e.g., of the order of the Alfven speed, the situation 
may change. For example, with Td yn — £/ f v a one gets / ~ Td yn /T pp ~ 10 independent of the 
proton energy (for B ~ 3 /xG, t ~ 3kpc, and n ism ~ 1 cm -3 ). In this case one would expect 
n p (E p ) ~ t pp x Q p (E p ) oc E p Sp , requiring thus a steeper injection index of s p ~ 2.75. 

Independent of which CR proton propagation model is the correct one, the production 
rate of the secondary pairs will depend on the observed spectrum of the CR protons (which 
we assume to be the same throughout the Galactic disk): 

Q e ±(E e ±) ~ r^ 1 n p {E p ) f e± {E e± /E p ) , (11) 

where f e ± (E e ± / E p ) is the num ber of pairs wi t h ene rgy E e ± produced by a CR proton of energy 



E p . Detailed calculations by iKelner et al.l (120061 ) show that for 1 TeV < E p < 1 PeV, the 



function f e ± (E e ± / E p ) is strongly peaked for E e ± / E p ~ 0.07 at the level / e T a ^ / e ± (0.07) ~ 4. 
As a result, the injection function of the secondary pairs should follow the energy spectrum of 
CR protons, namely Q e ±(E e ±) oc E e l v '. Since the secondary pairs obey the same transport 
equation as the primary electrons, their spectra can be calculated as discussed in §2. In 
particular, for high energies the escape term (as well as the Coulomb and bremsstrahlung 
energy losses) can be ignored giving n e ±(E e ± > lOGeV) — T rad x Q e ±(E e ±) ex Ej p for 
T ra d — Trad, t ■ More generally, the expected secondary pair to total electron ratio should 
vary with the energy roughly as 

r rad {E e ) AJ p {lAE e ) 



T pp Je\E e ) 
T^ism \ ( V'star \ ^ I E e 



pp 'PP 

0.75 



cm 3 / VeVcm 3 J V GeV 



(12) 



where J P (E P ) ~ 2.2 x 10 4 (E p /GeV)- 2 - 75 GeV" 1 m" 2 s" 1 sr" 1 is the observed CR proton flux, 
and J e (E e ) ~ 155 (i? e /GeV)~ 3 GeV -1 m~ 2 s _1 sr _1 is the observed CR electron flux. The last 
line in the above equation assumes we are in the Thomson regime with r ra d ~ r ra d, t oc E~± 
so that this ratio becomes n e ±/n e oc E~ ^ 5 with the particular value (n e ±/n e )i oGeV — 0.04 
for the assumed starlight density of u s t ar — 3eVcm~ 3 . This is in a disagreement with the 
PA MELA results indica ting; the fraction increasing with energy up to (n e ±/n e )i oGeV > 



0.1 (jAdriani et al.ll2009l ). However, as discussed above, at higher energies we are in the KN 

1/2 

regime, where r rad ~ T rad)KN oc E e , which will give rise to a flatter energy spectrum of 
the secondary pairs and hence to n e ±/n e oc E®± 5 . Note also that, since r rac i t T oc u~ tar , in 
regions of low (high) radiative field densities the expected e ± fraction will be higher (lower) 
for a given J p (E p ) and J e (E e ). 



13 



3.2. Photo-Pair Production 

One possibility for increasing the pair fraction in the CR spetrum is to introduce an 
additional, flatter spectral component consisting solely of the e* 1 pairs that outnumber the 
secondaries resulting from the proton-proton interactions^. However, this population cannot 
extend up to E e > 1 TeV energies, since this would viola te the high-energy cut- off measured 
in the CR electron spectrum by the HESS experiment (jAharonian et al.ll2008l ). A possible 
source of pairs that satisfy these requirements may be due to photon-photon annihila tion of 
TeV-energy 7-rays on starlight (jAharonian fc Atoyanlll99ll ; iMastichiadis et al.lll99ll ). The 
cross-section for this process has a sharp peak when photon energies satisfy the condition 
£o£ 7 = 2m^c 4 . Thus, the annihilation of Eq ~ e star ~ leV and e 1 ~ 0.5 TeV 7-ray pho- 
tons will inject into the ISM a relatively narrow energy distribution of pairs at the rate 
Q e ± i77 (i? e ±) oc 5(E e ± — m 2 e c A / e s tar) ■ Such a distribution cooling radiatively according to 
equation (j3J) will produce a flat-spectrum n e ±(E e ±) oc E~± (in the Thomson regime, or even 
a flatter one in the KN regime), instead of oc E~±' 75 expected for the secondaries resulting 
from the decay of ^ generated in the proton-proton interactions, as described above. 



However, the question is whether there will be sufficient number of such pairs to account 
for the PAMELA observations. In order to address this issue we use the ^-function approx- 
imation for the photon-photon annihila tion cross section a 77 (e ,e^) ~ (1/3) ut^o S[s — 
(2mgC 4 /e 7 )] (IZdziarski fc Lightmanlll985l ). from which we can calculate the absorption coef- 
ficient tt 77 (£ 7 ) — J m c 2 d,6o n (s ) cr 77 . If we also approximate the energy density of the soft 
(starlight) photon field by a monoenergetic distribution with total density n star and energy 
e star , namely n (e ) = n star 5(e - e star ) such that / de u (e ) = u star = n star e starj then the 
opacity becomes a 
optical depth to be 



77 \£ 1) 



<Jr/3) u star S[e s tar — ( 2m l c / £ j)]- From this we can evaluate the 



^ x 5 



2m?c 4 



E-star 



£0 



where 



77 3 



1 (^T Ustar £ star - 2 X 10 



-3 



Ustar 



3 kpc / VeV cm 



star 

~eV~ 



(13) 



(14) 



Since riL <C 1, most of the 7-rays freely escape the Galaxy and thus their number density 



77 



per energy is n, 
below. 



7V^7J 



(e/ c ) q 



where Q 7 (e 7 ) is the 7-ray production rate discussed 



3 Note that such a population cannot be accompanied by the additi onal population of proto n-antiproton 
pairs, since this would violate the observed proton-to-antiproton ratio ( Moskalenko et al. 20021 ). 
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3.2.1. Tertiary Pairs from Hadronic Interactions 

Proton-proton interactions, in addition to producing secondary pairs, also produce 7- 
rays (from 7r° decay) of similar spectrum and comparable intensity. These 7-rays could be 
the source of the pairs (which may be called tertiary pairs) in the above scenario. The 
rate of such 'hadronic' 7-ray production may be approximated as 

Qy,Pp(£~f) — Tpp n p(Ep) fy^Ey/Ep) , (15) 

where f^(£ 7 /E p ) is the number of photons with energy e 7 produced in a single proton-proton 
collision involving a C R proton with the e nergy E p . Just as in the case of secondary pair 
production, we refer to iKelner et al.l (120061 ). who showed that in the range 0.1TeV< E p < 



IPeV the function / 7 (e 7 /£ , p ) is peaked for E^/Ep ~ 0.1 at the level f™ ax ~ / 7 (0.1) ~ 6. 
The total production rate of such tertiary e ± pairs (with energies E e ^ £7/2) may 



be obtained from Q e ± ni (E e ) ~ 4c )| =2E (jCoppi fc Blandfordlll990l ). As before, 

inserting this in equation and carrying out the integration, we get the density ratio of 
tertiary pairs to total electrons at the same energy E e in terms of the observed CR proton 
and electron flux ratio, 



n P ± 



77/pp ^~PP E e E s tar Je(Ee) 



10" 



ri; 



3 kpc J Vcm 3 > 

where for the the bottom line we have assumed u star — u tot and used the Thomson regime 
for T rad {E e ). 

First, we note that because rfL 1 (see equation [T3l) the expected number of (tertiary) 
pairs from photo-pair process will be lower than that of the (secondary) pairs from proton- 
proton interaction. Second, because r° oc u star and in the Thomson regime r rad cx u~ tar , 
the ratio of photo-pairs to primary electron is independent of the energy density of the soft 
photon field, as long as it dominates over the other Galactic photon fields and the magnetic 
field. However, more importantly, this ratio increases with increasing starlight energy as 
n e ±/n e oc E^lal, and (in the Thomson regime) it increases linearly with electron energy. 
Therefore, in regions of the Galaxy containing high energy (ultraviolet) photons and for 
high energy electrons the photo-pair production may become important and even dominant 
(see below). Of course, the above result again will be modified by the inevitable KN effect. 
In this contex it should be emphasized that because of the flatter injection function of the 
tertiary pairs resulting from the photon-photon annihilation the KN effect should be more 
pronounced for them than for the primary electrons or the secondary pairs originating from 
the decay of due to proton-proton collisions. 
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3.2.2. Tertiary Pairs from Leptonic Interactions 

Yet another source of 7-rays which may annihilate on the starlight photon field and 
create additional e ± population is provided by the IC emission of CR electrons themselves. 
In order to estimate the expected relevance of this process, we need the rate of production 
of 7-rays, Q 7 ,j c ( £ 7 )> which will take the place of Q 7 , pp (£ 7 ) specified in the previous section. 
The IC rate is related to the IC emissivity ji C (£ 7 ) as Q 7) j C (£ 7 ) = 47rjj C (£: 7 )/e 7 , which can 
be obtained from the standard relation [s^ji C (e 7 )] ~ (1/4tt) [E 2 ic n e (E ejic )}/r rad (E e>ic ). Here 
E e> i C is the energy of electrons emitting 7-ray photons with energies e 7 , while r rad (E etic ) 
includes only IC cooling due to soft photons of energy e star . This gives 

Q,Ae,)=(^) 2 y%^- y (17) 

V 7 / rad\l- J e,ic) 

which replaces the photon production rate given above in equation (fl5l) . 

Following the same procedure as above we can evaluate the density of 7-rays, the rate of 
production of e ± pairs Q e ± tlJ (E e ±), and then the density of pairs in the ISM. In the Thomp- 
son regime e 7 ~ (4/3) E\ ic e star j m 1 ^ and thus (E £jic /e 7 ) 2 = 3/8 for e star = 2mlc A /e 7 , while 
in the KN regime E e ^ c /e.y ~ 1. As we will see below, for relevant CR energies we are closer 
to the KN regime so we will ignore the factor 3/8. We then obtain 



n, 



m 2 e c A [Je{.E e )/T rad )} Ee=2m 2 ci/es 



'yy/ic 



7 E e e 

star 



T 
> 



V'star \ ( E, 



x lO" 6 — ^ (18) 

V3kpcy VeVcm~ 3 / VGeVy V ; 

where the last line is evaluated for the Thomson regime with r rad oc E~ x . As evident, in this 
regime and for the observed J e (E e ) we obtain (n e ±/n e )i oGcV — 2 x 10~ 3 for £ ~ 3kpc and 
u star — 3eVcm~ 3 , independent of the soft photon energy. This implies that the production 
of TeV-energy 7-rays via IC emission of CR electrons — if proceeding in the Thomson regime 
- may dominate over the one resulting from the protons-proton interactions. On the other 
hand, the KN effect are expected to reduce the IC emissivity of ultrarelativistic e ± pairs 
within the consider photon energy range, and therefore both hadronic and leptonic processes 
may be in fact comparable. For the choice of model parameters appropriate for the average 
ISM conditions, this is not enough to account for the high positron fraction found in the CR 
spectrum. 
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3.3. Energy Spectra of Secondaries and Tertiaries 



As in case of primary electrons we now carry a more accurate determination of secondary 
and tertiary pairs by omitting most of the approximations used above. The production rate 
of the secondary pairs from proton-proton interactions is now obtained from 



1 



VP 



dE,, 



E 



- Jp(E p ) f e (E p , E e ) 



(19) 



where we use the analytic approximation for the function f e {E p , E e ) as given in lKelner et al.l 
J2006h . and fix J p (E p ) ~ 2.2 x 10 4 (E p /GeV)- 2 - 75 GeV" 1 m" 2 s" 1 sr -1 . For the production 
rates of (tertiary) pairs generated from annihilation of high-energy 7-rays with density n 7 (e 7 ) 
by the soft Galactic photon fields, we write analogously 



Q77 \E e 



4 

3 



-a T c u rad {e)\ £=mlci/E ^ n 7 ^ 7 J 



=2E ± ' 



(20) 



where again we have used the delta function approximation for the photon-photon anni- 
hilation cross section as before. For the expected small optical depth of photon-photon 
annihilation (r 77 <C 1, see above), the spectrum of 7-rays resulting from the proton-proton 
interactions is given by 



cr, 



VP 



dE p 
E n 



jpi.Ep) fj(E p , £ry) 



(21) 



with the function f 7 (E p ,Ej) denoting the number of photons with energy e 7 produced in 
a single collision involving ultrarelativistic proton wit h energy E p . Aga in, here we take 



the analytical approximation for fy(E„,e~) as given in iKelner et al.l (120061 ). noting that for 



E p ~ 0.1 TeV —1 PeV this may be further approximated by a simple function 

/ 7 (z) ~ 2.5 x- 1 exp [-9 a; ' 83 ] , 
with x = £j/E p (see in this context Hillasl [2005! ) . 



(22) 



Finally, for the case of 7-rays resulting from the IC emission of ultrarelativistic CR 
leptons, we calculate the appropriate photon energy spectrum as 



4vr£ 



Tly/i 



ic\<-"y J 



■J 



(23) 



where the IC emissivity jj C (e 7 ) is related to the observed electron flu x using the standard 



IC formulae with the KN effect included (IBlumenthal &: Gouldlll970l ). We fix this flux as 
J e (E e ) ~ 155 (E e /GeV)- 3 GeV" 1 m~ 2 s~ x sr" 1 (and cutting-off exponentially at E e = 2 TeV). 
Inserting then the resulting e ± pair production rates in equation (j3J), we obtain the individual 
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and total e ± pair fluxes Jjg(E e ) = J p JL(E e 



>+j]2 /pp 



discussed above and for the same seven parameters s e , 
used in calculation of the primary electron spectra. 



(E e ) for all three mechanisms 
£ and B 



Ee : maxi V'dusti ^stari H". 



FigureH] shows the energy spectra of secondary leptons produced in proton-proton col- 
lisions, J p ±(E e ) (solid lines), and of tertiary pairs produced via absorption of high-energy 
7-rays generated in either hadronic or leptonic processes (dashed and dotted lines, respec- 



tively), r p 2 /pP (E 



,,) and J]2^ lC {E e ), for the same set of parameters used in Figured] with same 



colors, except for the magnetic field set at B = 1 /j,G. As shown, the cases with a weak 
starlight but strong dust emission are quantitatively similar to the cases when the dust 
and starlight energy densities are comparable. Only in the cases when the ratio u sta r/udust 
is high, the KN effect flattens the energy distribution of secondary leptons resulting from 
proton-proton interactions significantly, and only at high (E e > 100 GeV) energies. In all 
the above cases, however, the direct e ± pair production in the proton-proton collisions is 
the dominant source of the positrons, and the contribution of the other two processes (i.e., 
of the tertiary pairs) to the positron flux is less than 1% except at high energies where it 
could reach 10%. The contribution of secondary and tertiary electrons to the observed elec- 
tron spectrum is even less, being on the order of < 10% and < 0.1% for proton-proton and 
photo-pair processes, respectively. 

The above result is illustrated in Figure[5j where we compare with different measure- 
ments the computed ratio of (both secondary and tertiary) positron and electron fluxes, 
normally denoted in the literature as 0(e + ) and (f>(e~), respectively, 



+ 



2 [j e (E e )/Jg(E t 



(24) 

3eVcm" 3 , n ism = 
= 2.75TeV). Here 

red symbols denote the PAMELA data (lAdriani etal.l 120091 ). cyan and blue symbols the 



for the same model parameters as considered in Figured (namely u s t ar 



1 cm 



u dust = O.leVcm" 3 , B = 3/xG, I = 3kpc, s P = 2.42 and E, 



Boezio et al 



( I Golden et al. 



HEAT data (Barwick et al.l ll997l: iBeatty et all 120041 ). yellow symbols the CAPRICE data 



e, max 



2000), and the magenta ones the measurements with the imaging calorimeter 
1994! ) . As evident, due to the KN effects and inclusion of tertiary pairs the 
^ fraction decreases only by a factor of 2 between E e ~ 10 Ge V and 200 GeV. Even thoug h 



this is a much less rapid decrease than typically expected (see Moskalenko fc Strong! Il998l ). 
the PAMELA results in the high energy [E e > 20 GeV) range cannot be reproduced with 
our conservative choice of model parameters. 
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3.4. High Positron Fraction 

The above results show that it would be rather difficult to increase the fraction of sec- 
ondary pairs just by the photo-pair processes for the average ISM conditions. One needs 
different conditions for the production of the relatively high positron-to-electron fraction in 
the CR spectrum which also increases with energy, as claimed by the PAMELA experiment. 
In our model, those would require increasing the energy density of the starlight emission up 
to u s tar ~ 300eVcm~ 3 , and of the ISM number density up to rii sm ~ 80 cm -3 , keeping at 
the same time relatively low level of dust emission (udust ~ O.leVcm -3 ) and magnetic field 
strength (B < 10 //G). FigureE] (bottom panel) shows the e ± fraction expected for such a 
choice of model parameters, which agrees with the PAMELA data within the energy range 
not affected by the solar modulation. The corresponding total electron spectrum is com- 
pared with the Fermi and HESS data in the top panel of Figured As evident we get again 
a very good agreement but now we need an even steeper injection spectrum of the primary 
electrons (s e ~ 2.65). 

The set of model parameters considered in Figure[6] should be regarded as illustrative one 
only, not necessarily being justified for the local ISM. We note, however, that it corresponds 
to the optical depth for annihilation of Galactic 7-rays formally less than (though close to) 
unity (see equation[TJJ , and to a small (< 0.1) ratio of number densities of 7-ray photons 
and CR electrons with the same energy e, as required. In fact, we have 

n 7 (e) _ n JtPP (e) _ Q£ J p (10e) c _ in _ 4 ( n ism \ ( e \ - 25 



n e {e) n e (e) cr pp J e (e) \cm 6 / VGeV/ 

(see the discussion in section 3.2 above). 

The invoked increased level of the starlight energy density and of the gas number density 
could be more appropriate around supernova remnants where the injection of the Galactic 
CRs is taking place. Hence, the results of our analysis may indicate that ultrarelativistic 
particles generated in the Galaxy undergo most of their interactions near their sources, but 
propagate m uch more freely from these regions to th e Earth (see in this context recent 



discussion in lHigdon et al.ll2009l ; ICowsik fc Burchl 120091 ). In fact it may be sufficient if high- 



energy positrons, but not necessarily electrons, are trapped in the regions characterized by 
the enhanced photon and gas densities. There may be even physical justification for such 
a situation. For example, CR protons streaming along large-scale magnetic field in the far 
upstream of supernova shocks with super- Alfvenic speed ma y excite resonant Alfven waves in 



a form of coherent cir cularly-polarized cyclotron radiation (jLerchd 119671 ; iKulsrud fc Pearce 



19691 ; ICesarskyi I198CH ) . Due to the particular helicity of the generated waves, they will 
interact with positrons of gyroradii comparable to their wavelenghts (i.e., to gyroradii of CR 
protons generating the turbulence), but not with the electrons. As a result, the electrons 
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will propagate much more freely along the Galactic magnetic field to the Earth, experiencing 
the 'average' ISM conditions^. The results presented in Figure[3] regarding the observed CR 
electron spectrum would then be appropriate. CR positrons, on the other hand, will undergo 
enhanced scattering in vicinities of their sources resulting in their increased fraction in the 
observed CR spectrum around 100 GeV energies, as presented in Figure[6] (bottom panel). A 
quantitative description of such a possibility would require different treatment of the positron 
and electron transport within the Galaxy. This is beyond the scope of this paper. The point 
is, however, that the efficient trapping of TeV-energy CRs in vicinities of supernova remnants, 
either charge- dependent or not, may justify the high values for the starlight and gas densities 
invoked to explain the PAMELA data in a framework of our model. 

Let us mention in this context that in the local environments of SNRs additional pro- 
cesses may operate leading to an increase in the CR positron-to-electron fraction. These 
include enhanced interactions of freshly accelerated CR protons with an intense high-energy 
photon field of young remna nts, generating t hus additional secondary pairs via the photo- 
mezon production process (IHu et all 120091 ). or the direct a cceleration of secondary pairs 
injected into the immediacy of SNR shocks via pp collisions (]Blasill2009l ). 



Yet one more process which may be relevant in the discussed context is the creation of 
; pairs by photons in the electromagnetic field of ultrarelativistic electrons, referred in the 



Mastichiadis et al. 


1986; 


Mastichiadis 


1991 



Dermer fc Schlickeiserl Il99ll ). This process occurs when the energy of the incident photon 
in the electron rest frame exceeds 4 times the rest energy of the electron, e' > e cr = 4m e c 2 . 
For the starlight parameters considered in this paper, namely e ~ leV, this criterium is 
marginally fulfilled only in the 'head-on' interactions with the highest energy electrons, E e ~ 
ITeV, since only in such a case e' ~ 2 e E e /m e c 2 ~ 2e cr . For all the other angles between 
the direction of an interacting electron and starlight photon, and for all the lower-energy 
electrons, we have obviously e' <C e cr . Nevertheless, the TPP may be of a primary importance 
if the soft photon energies are higher than anticipated here, say e ~ 10 eV. Then the head- 
on collisions of such UV photons with the TeV-energy elec trons will produce effectivel y 
e ± pairs with energies E e ± ~ 0.5 (E e /e 



\ l l 2 m e c 2 ~ 0.1 TeV (see lDermer &: Schlickeiserl Il99ll ) 



i.e. exactly within the energy range of the PAMELA excess. Note that the energy losses 
of thus produced pairs should be dominated by the IC scattering deep in the KN regime, 
and hence the spectral pile-ups discussed in this paper will flatten additionally the injected 
positron spectrum around E e ~ 10 — 100 GeV energies. That is because the TPP cross- 
section, otpp ~ Oif s x or, exceeds the IC cross section (due to the KN supression of the 



4 Note that the returning current will be assured by the ambient plasma, and would involve sub-thermal 
bulk velocities of ISM particles due to the expected high number density of ISM within the Galactic disk. 
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latter one) only for e' > 300 m P c 2 , while the TPP cooling rate exceeds the IC cooling rate 
only for e' > 10 5 m e c 3 (jMastichiadia Il99ll ; iDermer &: Schlickeiserl Il99ll ). As a result, if the 
sources of Galactic CRs are associated with an intense UV photon field, the most recent 
PAMELA results may be possibly explained with much less extreme ISM parameters than 
discussed in this section. 



4. Summary and Discussion 

In this paper we show that the observed excesses in the energy distribution of the 
Galactic CR electrons around energies E e ~ 0.1 — 1 TeV may be easily re-produced without 
invoking any unusual source of ultrarelativistic electrons (or e pairs), such as dark matter 
annihilation/ decay or some nearby astrophysical object (e.g. a pulsar), other than the general 
diffuse Galactic components of CR electrons and protons injected by supernova remnants. 
The model presented here assumes an injected spectrum of electrons (power-law with index 
s e ) and evaluates their observed energy distribution based on a simple and most commonly 
invoked kinetic equation describing the propagation of CR electrons in the ISM. The main 
process affecting this outcome is the cooling of the injected electrons by their interaction 
with the ISM photons (via IC scattering). The interactions of electrons with ISM turbulence 
produces negligible re-acceleration and determines their escape time. The escape timescale 
also turns out to be somewhat longer than the cooling time in the relevant range of electron 
energies. The new physical effect that is the source of the observed excess is the Klein- Nishina 
suppression of the IC cooling rate, which becomes important right around TeV energies. 
With a very reasonable choice of the model parameters characterizing the local interstellar 
medium (u star ~ 3eVcm -3 , Ud us t ~ u cm b ~ 0.3eVcm~ 3 , B ~ 3/xG, and n ism ~ lcm -3 ) we 
can reproduce the most recent, and perhaps the most reliable observations by Fermi and 
HESS, but not the sharp feature claimed by ATIC. Interestingly, in our model the injection 
spectral index of CR electrons becomes comparable to, or perhaps equal to that of CR 
protons, namely s e ~ s p ~ 2.4. 

The Klein-Nishina effect will also affect the propagation of the secondary pairs and 
can produce deviations from a power-law in the observed spectra of such pairs. In particular, 
it can affect the positron-to-electron ratio. We have explored this possibility by considering 
two mechanisms for production of e pairs (and therefore positrons). The first is production 
of pairs due to the decay of tt^s generated by interaction of CR nuclei with ambient protons. 
The second source discussed here is the pair production due to annihilation of diffuse Galactic 
7-rays interacting with the starlight photon field. We consider two sources of the Galactic 
7-rays. The first is related to the decay of 7r°'s also produced in proton-proton interactions 
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and the second is due to the IC scattering of primary CR electrons by the diffuse Galactic 
photon fields. We show that indeed there will be deviations from a simple power-law in the 
spectra of thus created pairs (as well as in the positron-to-electron flux ratio), similar 
to the observed one. However, the relatively high observed positron fraction that increases 
quite steeply with energy, as observed by PAMELA, cannot be explained by the conservative 
set of the model parameters used above, which corresponds to the average values expected 
in the Galactic disk. We can however reproduce the PAMELA result by increasing the 
energy density of the starlight photon field and of the ISM number density up to the levels 
Ustar ~ 300eVcm -3 and rii sm ~ 80 cm -3 . With these new values we can also fit the Fermi 
and HESS data, though with somewhat steeper injected spectrum of the primary electrons 
than required before (s e ~ 2.65). 

The required increased level of the starlight energy density and of the gas number density 
may be regarded as unlikely for the local interstellar medium. However, such a choice of the 
model parameters could be more appropriate around supernova remnants where the injection 
of the Galactic CRs is taking place. A possible solution to this problem may be that CRs 
undergo most of their interactions near their sources, being efficiently trapped thereby by self- 
generated CR-driven turbulence. Interestingly, such a trapping may be charge- dependent, 
affecting positrons more than the electrons. A possible cause of this could be if the domi- 
nant CRs, namely protons, generate Alfven waves of a particular helicity which scatter and 
therefore trap positrons more efficiently than electrons in the regions characterized by the 
enhanced photon and gas densities. Alternatively, higher than considered here energies of 
photons associated with CR sources may reduce significantly the invoked 'extreme' values 
of the model parameters, due to even more severe KN effects and additional (triplet) pair 
production processes expected to occur in an intense UV radiation field. 

We note in this context that the qualitatively similar effects to the ones analyzed here for 
the case of our Galaxy have be en discussed pr e yiousl y for the case of the host galaxy of nearby 



radio source Centaurus A by IStawarz et al.l (120061 ). The theoretically predicted isotropic 



galactic-scale halo of ultrarelativistic pairs thereby (with the energy distribution shaped 
by the KN and 7-ray annihilation processes), and in particul ar the resulting TeV em ission, 



has been possibly already detected by the HESS instrument (lAharonian et al.ll2009b[ ). 
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Fig. 1. — Left panels: Models of the target photon fields including the CMB and different 
valuse of starlight and dust emission. Each panel has three curves for Udust and/or u sta r equal 
to 0.3, 1, or 3eVcm~ 3 , varying independently or together (black, red, and blue curves). In 
the top panel u star is set as 0.3eVcm~ 3 . In the middle panel Udust is set as 0.3eVcm -3 . 
In the bottom panel u sta r = Udust- Two different values of the magnetic field densities are 
also shown: B = 1 fiG (solid horizontal lines), and 3 fiG (dashed horizontal lines). Right 
panels: The energy dependence of the energy losses timescales (multiplied by energy) for 
the Galactic CR electrons corresponding to the different levels of the Galactic photon and 
magnetic fields shown on the left panels, and n ism = 1 cm -3 . 
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Fig. 2. — The energy spectra of primary electrons corresponding to two different injection 
spectral indices, s e = 2.0 (left panels) and 2.2 (right panels), for the same set of the 
model parameters given in Figured] 
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Fig. 3. — Comparison of the observed spectra of the Galactic CR electrons with model 
spectra calculated for u star = 3eVcm~ 3 , n ism = lcm~ 3 , Ud us t = O.leVcm -3 , B = 3/zG, 
£ = 3kpc, s e = 2.42 and E e ^ max = 2.75 TeV. The black solid line corresponds to the energy 
spectrum of the primary CR electrons calculated usi ng equations (4-7) . The data points 
corr espond to differen t measurements by ATIC (black; 



low; Torii et al. 



2008) , emulsion chamb ers (magenta; 



Chang et alJl2008h. PP B-BETS (yel- 



Kobavashi etal]l2004h. HE SS fblue 



and cyan; lAharonian et al.ll2008l . l2009al . respectively), and Fermi (red: lAbdo et al.ll2009l ). 
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Fig. 4. — The energy spectra of secondary leptons produced in proton-proton collisions (solid 
lines), and of tertiary pairs produced via absorption of high-energy 7-rays generated in either 
hadronic or leptonic processes (dashed and dotted lines, respectively), for the same set of the 
model free parameters as given in Figured] (black, red, and blue curves on different panels), 
except for the single value of the magnetic field B = 1 /zG. 



-30- 



0.300 



0.200 

" 0.150 



^ 0.100 



+ 



0.070 



^0.050 



0.030 



0.020 




5 10 20 
E e [GeV] 



100 200 



Fig. 5. — The positron-t o-electron ratio from different me asurements by PAMELA (red sym- 



bols; lAdriani et al.ll2009l ). HEAT (c yan and blue symb ols; iBarwick et al.lll997l ; iBeatty et al. 



20041 ). C APRICE (yellow sy mbols; iBoezio et al.ll2000l ). and imaging calorimeter (magenta 



symbols; iGolden et al.lll994l ). compared with the model result (line) for the same model 
parameters as considered in Figure[3j namely u sta r 
O.leVcm" 3 , B = 3fiG, t = 3kpc, s e = 2.42 and E e , 



3eVcm 3 , n ism 
= 2.75 TeV. 
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Fig. 6. — Top panel: The total energy spectrum of cosmic ray leptons calculated for u star 



300eVcm 



-3 



80 cm 3 , Udust = O.leVcm 3 , B = 3/uG, t = 3kpc, s 



e, in] 



2.65 and 



E e ,max = 1.55 TeV (assuming super-exponential cut-off in the electron injection function). 
The data points correspond to the measurements by HESS (blue and cyan symbols) and 
Fermi (red symbols). Bottom panel: The resulting positron-to-electron ratio compared 
with the measurements by PAMELA (red symbols). 



